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Abstract 

We demonstrate theory and computations for finite-energy line defect solutions in 
an improvement of Ericksen-Leslie liquid crystal theory. Planar director fields are 
considered in two and three space dimensions, and we demonstrate straight as well as 
loop disclination solutions. The possibility of static balance of forces in the presence of 
a disclination and in the absence of flow and body forces is discussed. The work exploits 
an implicit conceptual connection between the Weingarten-Volterra characterization of 
possible jumps in certain potential fields and the Stokes-Helmholtz resolution of vector 
fields. The theoretical basis of our work is compared and contrasted with the theory 
of Volterra disclinations in elasticity. Physical reasoning precluding a gauge-invariant 
structure for the model is also presented. 

1 Introduction 

While liquid crystal theory, especially that of nematics, is by now a well-established branch 
of condensed matter physics |KL03 j. there is one aspect in which the classical theory may 
be considered deficient. Disclination line defects are an integral part of the physics of liquid 
crystals - however, there does not exist a classical theory, whether Oseen- Frank for statics 
or Ericksen-Leslie in dynamics (which naturally subsumes the static theory), that predicts 
bounded energy in a finite body containing a ±|-strength disclination line. It is perhaps for 
this reason that it is commonplace amongst workers in liquid crystal theory to associate a 
line defect only with a singularity, e.g. the 'escape' solution for a +1 strength line defect 
jCK72j is often not referred to as a defect solution. 

Director fields with infinite energy are practically problematic from the point of view of 
finding solutions - e.g. in statics, energy cannot be minimized in a class of functions where 
all competing fields with defects have undefined (infinite) energies. Practical difficulties 
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remain in dynamics, more so if the goal is to specify dynamics of defect lines. While for the 
purpose of static Oseen-Frank theory involving a few defects one overcomes such problems 
by excluding from analysis a core region of small volume along the defect line where one 
assumes the theory does not hold and ascribing a finite energy to it from some unspecified- 
in-theory nonlinear effects, for situations where large numbers of defects may be involved, 
such an approach is untenable. To quote Ericksen [ Eri91j "I am interested in seeing the 
development of a mathematically sound theory of defects which might be at rest, or moving. 
...serious difficulties are encountered with some observed kinds of disclinations (line defects). 
According to the aforementioned theories [Oseen-Frank, Ericksen-Leslie] , these more violent 
singularities cause energy integrals, etc., to diverge. In dealing with rather specific situations, 
workers have patched up the theory by excluding a tube of small radius, assigning a finite 
"core energy" to it. I just do not believe that one can use such ideas as a basis for developing 
satisfactory mathematical theory, particularly for moving defects. Such phenomena are of 
interest for both nematic and cholesteric liquid crystals. In the latter, they are associated 

with the interesting "blue phases" " 

In jAD12j a theory for the dynamics of non-singular line defects has been proposed as 
an extension of Ericksen-Leslie (Les92j theory. The model works with augmented kinematics 
involving the director field and an incompatible director distortion field as a replacement for 
the director gradient field. The main point of departure is that the director distortion is 
not curl-iiee^ in general. The curl of the field has the kinematic meaning of an areal den- 
sity of lines carrying a vectorial attribute and is referred to as the director incompatibility 
field. For the 1-constant Oseen-Frank energy approximation, it was shown in |AD12| that 
the developed formalism can indeed predict the outside-core, energy density distribution of 
Frank's planar director distributions for line defects with strength given by all integer mul- 
tiples of |. Moreover, all such distributions have finite energy in the whole body (including 
core). The core is represented by precisely the region where the director incompatibility 
field is non- vanishing. However, it was not shown in that work how to extract a physically 
realistic continuous director field resembling Frank's planar defect fields on the whole body, 
and neither was it shown as to what form the director distortion tensor might take within 
the theory to represent planar distortion fields of disclination line defects. In this paper we 
achieve that goal. 

The modeling of bounded energy line defects in statics started with the pioneering work 
of jCK72j within classical Oseen-Frank theory. The same question served as motivation 
for the development of the variable degree of orientation extension of classical Ericksen 
Leslie theory in [Eri91j . Within the De-Gennes Q tensor formalism it has been pursued in 
jSKH95j . jKVicv99] . |BZ07j . jMZlOj . In the dynamic case, representative studies are those 
of jHKGLGCCOS) and |YFMW09j within the Q-tensor approach and [LLOOl ILWOOj within 
the Ericksen-Leslie approach. In |LLOO[ ILWOOj non-singular line-defect solutions appear due 
to relaxation of the unit-magnitude constraint on the director field. In contrast, in our work 
the demonstrated finite energy defect solutions satisfy the strict unit-vector constraint on 
the director field. 

The defect solutions we develop here are finite energy analogs of classically accepted 
'energy-minimizing' (infinite total energy) singular solutions of Oseen-Zocher-Frank theory. 
As is customary in the liquid crystal literature when discussing statics of defects solutions, 
we (almost) ignore balance of forces. What we are able to show in this regard (Section [5]) is 
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that static balance of linear momentum is satisfied outside the core(s) of disclination(s) in 
the absence of a body force within our model (assuming an incompressible fluid). Moreover, 
when there is a single disclination in the body, the natural boundary conditions for static 
force and moment balance representing no applied forces and couples on the boundary imply 
that the resultant force on any surface enveloping the disclination core vanishes. Within the 
core, we discuss possibilities afforded by the |AD12j theory in satisfying balance of linear 
momentum without external body forces or flow. 

We consider this work to be primarily a characterization of the kinematic augmentation 
of the Ericksen-Leslie theory achieved by the [AD12J model, since results are not presented 
regarding the connection of the presented solutions to (quasi) equilibria (cf. |CP89j ) of the 
latter. 

This paper is organized as follows: in Section [2] we provide notation. Section [3] describes 
the analytical basis of our work; Section [4] describes the flnite-element based numerical 
scheme employed to compute non-singular defect solutions and the results. In Section [5] 
we discuss static balance of forces for the developed solutions, and Section |6] ends with 
some concluding remarks. Two appendices provide explicit solutions that are necessary for 
discussing our results. 



2 Notation 

A represents the inner product of two vectors, while represents the trace inner product 
of two second-order tensors, A : B = AijBij. 

In rectangular Cartesian coordinates and corresponding bases and components, for a 
vector fleld A and a scalar fleld 9 we write: 

div A = Ai^i 
{cutI A^i = Cijj^Aj^j 
{grad9)i = 

For a second-order tensor fleld A, we deflne 

(curl A). ^ = ejnjkAikj- 
The following list describes some of the mathematical symbols we use in this paper: 
n : director 

E : director distortion tensor 

: director distortion vector 
/3 : director incompatibility tensor 

: director incompatibility vector 
K : defect strength 
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: normal vector 


T 


: Cauchy stress tensor 


A 


: couple stress tensor 


^ 


: free energy per unit mass 


p : 


pressure 


e : 


angle of the director field 


A 


: layer field 


Ax 


: incompatible part of A 


grad z : compatible part of A 
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Theory 



We consider the static (i.e. director and positional inertia-less) governing equations of the 
framework introduced in |AD12j and demonstrate that there exists finite-energy states in 
the model that are solutions corresponding to line defects with planar director distributions 
within the Oseen-Frank constitutive assumption. 

3.1 Governing equations 

Static Balance of Angular Momentum is given by the statement 

Ericksen's identity for the theory (i.e. a necessary condition for frame-indiflFerence of the 
energy function) is the statement 

/dip ^^7? \ 

[d^''' + dE~ ~ '^'dE~ J ^ °' 

Static Balance of Linear Momentum is the statement 

T = 

For a free-energy density of the form ijj (n, the constitutive statements for the stress 
and couple stress tensors for the theory jAD12 j read as 

_ A - 
4 



A. 



T 



where p is the pressure arising from the constraint of incomressibihty. Using these statements 
along with the Ericksen identity in angular momentum balance with some manipulations 
yields the fundamental governing partial differential equation of our work: 



dip ^ , / 



0. (1) 



Note the interesting fact that when E^j = n^j, i.e. there is no director incompatibihty and 
hence no defects, we have the classical statement of angular momentum balance in statics 
|Ste04j . Inasmuch, this term in our theory may be interpreted as giving one explicit form to 
Ericksen's "internal body moments" (EriQlj . gf^, arising from the presence of defects in the 
body. 

We record the statement of balance of linear momentum, 

-p^,-(Er^^] =0, (2) 



to be used later in Section [H 

For the sake of this paper, we consider an ansatz consisting of director distributions 
parametrizable by a single angle field 9 on the body and director distortion fields of the form 

^:=^®< (3) 

where is a vector field with possibly non- vanishing curl. When there are no line defects, 
we require = gradO so that E = gradn. 

Let the core region for a given E^ field on the body be the set of points on which 

f3^ := curlE^ ^ 

and 

/3 := curl{E — gradn) ^ 0. 

We refer to an E^ field as containing an isolated line defect if its (3^ field is non zero only in 
a cylindrical region that forms a closed loop or extends from one boundary of the body to 
another. 

The goal now is to construct pairs of (£", n) fields 

1. that satisfy ([T]); 

2. that produce Frank's isolated line defect-like energy density fields outside the core 
region, while producing bounded energy in finite bodies (including the core region), 
when no moments are applied on the boundary of the body; 

3. that produce Frank's wedge disclination-like director distributions except the core and 
possibly another region of 'small', but non-zero, volume in the body; 

4. for which the line integral of iJ^, along any closed contour surrounding the core of an 
isolated line defect of strength X, evaluates to 2tiK. 
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3.2 Construction of bounded energy distortion and director solu- 
tions 

We work with a rectangular Cartesian coordinate system with unit vectors e^,z = 1,2,3. 
6 represents the angle of the director measured counter-clockwise from the Xi-axis in the 
Xi — X2 plane (looking down the X3-axis). Thus 

n = cos 9ei + sin 6^62- (4) 

For the sake of illustrating the essentials of our approach, we work with the 1-constant 
Oseen- Frank energy density approximation 

^ (n, E) = ]^kE:E = ^kE^.E,, = ^kE^E^ (5) 



and 
Noting 

and ([5])-([6]), ([T]) reduces to 



^Es, = ^^El (6) 



Erj - Tlrj = {Ej - 6j) 

dn 



= 0. 



Considering (4|) and realizing that eirsTir^Q^Oj = 0, static angular momentum balance re- 
duces to 

e. n ^E'.-O 

which is satisfied if and only if 

divE^ = 0. (7) 

To fulfill the fourth requirement on the field ^ let us now consider a distribution of 
f3^ := curlE^ in a right-cylinder parallel to X3 such that 



/ • 63 dxidx2 = 27rK, 
J A 



(8) 



where A is the set of points representing the cross section of the core cylinder. Thus, in 
addition to ([g]), we require 

curlE^ = f3^. (9) 

Finally, to impose the fact that we seek non-trivial solutions for an unloaded body, we impose 
vanishing moments on the boundary of the body which, under the current ansatz is satisfied 
if 

E^ -1^ = (10) 

on the boundary of the body with normal field 

It is straightforward to see that, for a specified field on a simply-connected domain, 
Q-g-lllOl) has a unique solution. In |AD12) . these equations are solved explicitly for a 
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radially symmetric /3 distribution (see Appendix |B|) , and the outside-core result shown to 
be exactly the same as the gradient of the Frank angle field for straight wedge disclinations. 
In the general case, (i.e. not necessarily radially-symmetric (3^ field, but still localized in core 
cylinder distributions) we would now like to extract the outside-core result into a gradient 
of a continuous everywhere scalar angle field 9 which will serve to define the n field of the 
pair (£", n) describing an isolated defect. Of course, in the process we would like to have 



stiU be the solution to 0-(|§-(|lO|) = grade when f3^ = 0. 

To perform this extraction, suppose for a moment it is possible to construct a square- 
integrable vector field A such that 

curlX = -f3^ (11) 



on the body, where the /3 field is identical to the one occurring in ([9|). Then a Stokes- 
Helmholtz resolution of the field A can be written as 

A = A^ + gradz 

with 

curl X± = curl A = -^^ (12) 
div A^ = 

and 

X± • 1/ = on boundary of body. 



and because of the uniqueness of solutions to ([7|)-([9|)-(10), = — A^. 

Suppose further now that the field A is non-vanishing only on a layer-like region of non- 
zero volum^thdit may be visualized as a terminating fattened 2-d surface which, moreover, 
contains as a proper subset the core cylinder on which has support (Figure [l]). 

Then, as an immediate consequence we have that 

—gradz = A^ = —E^ outside the layer. 

Thus, if we were to now declare the potential z as the required angle-field 0, we would have 
completed the needed extraction. 

In practice (Section [3]), we implement the above idea by defining 

E^ := grade -X (13) 

and requiring this combination to satisfy ([7])-([9])-(10) so that 



grade -X = -X^ (14) 

is enforced by uniqueness implying grade = gradz. When accompanied by the requirement 
that A = if = 0, the strategy ensures that E = gradn in the absence of defects in the 
body. We note the important fact that the Stokes-Helmholtz decomposition of an L2 vector 



^ Cf. |DeW73^ ; DeWit considers a distributional layer field in the context of disclinations in solids related 
to defects in positional order. 
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layer 




Figure 1: Schematic of cross-section of a 3-d body showing layer and core geometry. 



field implies that grad z is square-integrable and therefore grad 6 cannot have a non-square- 
integrable singularity in the body. The set of Definitions and (13) constitutes a primary 
result of this paper. 

So, the only question that remains is whether such a layer- field A satisfying (11) can in 
fact be constructed. As we show in Section |4| it suffices to consider the layer region as in 
Figure [l| prescribe an appropriate constant vector field (representative of an interpolation 
of the required jump in the value of 0, i.e. 27r/C, in the direction transverse to the layer) for 
the values of A outside the core in it, and tapering this interpolated jump to zero over the 
width of the core. Outside the layer, A is assigned to vanish. Since A represents a transverse 
'gradient' in the layer whereas in the core this transverse variation has a gradient in an in- 
plane direction, it cannot have vanishing curl in the core. Thus we represent a non-singular, 
but localized, core of a line defect. 

We now ask the question of the region where /3 is non-vanishing corresponding to an 
isolated wedge disclination. The existence of a differentiable extracted 9 field (which can be 
arranged as can be seen from the construction in Appendix [A]) implies that the /3 field corre- 
sponding to an isolated defect would be non- vanishing at most in the layer, since grad 9 = 
outside the layer and 

^2 77/ Qrii 

Our use of a layer field may be interpreted as a regularized analog of the possible termi- 
nating discontinuity in the displacement field over a 'cut' surface in the Weingarten-Volterra 
(WV) process of elastic dislocation theory, adapted to a much simpler situation than elastic- 
ity; thus, the construction shows direct links between the Stokes-Helmholtz (SH) resolution 
of a vector field (in this instance, the layer field) and the Weingarten-Volterra process. In 
particular, how the ^curV part of SH encodes, through a smooth field, topological content 
that can only be represented by a nasty singularity in the gradient of a discontinuous poten- 
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tial field in the WV process. And, how the 'gradient' part of SH of an appropriately designed 
vector field (the layer field) can represent most of the characteristics of the discontinuous 
potential field of the WV process away from the discontinuity. 

3.3 Does topological defect density determine the director distri- 
bution? 

The construction above lays bare an interesting fact. Note that is uniquely determined 
by the defect density field (and balance of moments) but gradO^ and hence the predicted 
director distribution, is defined by the compatible part gradz of A. An alternative way of 
seeing this is to substitute 

= gradO — — gradz 

into ([7|-([9l)-([l0l) (this idea is not restricted to the 1-constant energy in any way). Thus, while 
the energy density distribution, and therefore the couple stress, can be correctly predicted 
purely from the knowledge of the defect distribution (in this static setting), uniquely pre- 
dicting the director distribution requires additional physical input in the presence of defects. 
In particular, the layer-like A field consistent with a specified isolated defect density field 
orients the nonsingular director distribution of the defect, as is shown in Section [4} To draw 
an association with classical potential theory, specifying the surface of discontinuity of a 
scalar potential field whose gradient is required to match, except on the discontinuity sur- 
face, a prescribed, smooth, irrotational vector field on a doubly-connected region containing 
a toroidal or through-hole eliminates the vast non-uniqueness associated with such definition 
of a potential otherwise. The argument also reveals why the classical theory with only the 
director field as a degree of freedom is necessarily limited in the context of modeling defects. 

It is also interesting to observe that this situation is entirely analogous to the extension 
that was required for making a prediction of permanent deformation [AchOlJ in the elastic 
theory of continuously distributed dislocations [Kr oSlj . |Wil67j . Comparing and contrasting 
with prevalent notions in gauge theories, 6 and the total displacement field of plasticity 
theory are not inconsequential fields to be gauge-transformed away even when they have no 
energetic cost through the Lagrangian in the presence of gauge fields (A, in this context), 
but direct physical observables that participate in understanding energy dissipation, shape 
changes, and optical response of the body. Hence, their evolution requires additional physical 
specification beyond the topological defect density field. 

3.4 Relationship with disclinations in elasticity theory 

Our construction of the A field draws direct motivation from a standard question of potential 
theory related to understanding the allowable jumps of a scalar potential across a 2-d surface 
in a multiply-connected domain described by a through-hole in a 3-d body. For example, 
consider Figure [T] with the layer shrunk to zero thickness and also the consequent core to a 
line of puncture through the plane of the paper; the core-line represents the through-hole. 
The 2-d surface is such that a cut along it renders simply-connected the body with the hole. 
Additionally, the gradient of the potential field is required to match a prescribed irrotational 
vector field on the simply-connected domain obtained through the cut. 
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We note here that the question above is much simpler in its details than the question 
posed by Weingarten and Volterra related to understanding the jumps of a displacement 
(vector) field, when its symmetrized gradient is required to match a prescribed symmetric 
tensor field which satisfies the St.-Venant compatibility condition of linear elasticity theory in 
the simply-connected domain resulting from a planar cut of the multiply-connected domain 
with the hole. Indeed, in this latter case, a key construction is a skew-symmetric tensor field 
(infinitesimal rotation) from the prescribed symmetric tensor field, whose jump across the 
cut-surfaces (that serve as parts of the boundary of the simply-connected domain) relates 
directly to the disclinations of linear elasticity theory. It should be clear, then, that in the 
liquid crystal theory context, even if one chose to discuss the issue of allowable jumps in the 
director field context rather than the angle parametrization field, there is no natural way 
a relevant symmetric tensor field arises from the director gradient whose uniqu^ infinites- 
imal rotation field in a simply-connected domain may then be used to define the notion of 
an isolated singular disclination in liquid crystal theory. Given that the Oseen-Frank and 
Ericksen-Leslie theories are geometrically exact, the situation is not very diflFerent if a link 
is attempted between liquid crystal disclinations and the singular defects that characterize 
the jumps in finite rotation associated with an elastic strain field that satisfies the finite 
strain compatibility conditions in a multiply-connected domain. Again, the issue crucially 
hinges on the question being posed in terms of the basic tensor field being a positive-definite 
symmetric second order tensor field (the Right-Cauchy-Green tensor field) and there does 
not seem to be any natural way in which the director gradient may be considered to generate 
such a field. Such an observation is consistent with that made in |KF08j (Sec. I.E.I, p. 67): 
"This explains the interest in reconsidering the continuous theory of defects, although new 
concepts have to emerge. The case of mesomorphic phases requires an extension of the theory 
of continuous defects for solids to situations where there is locally only one physical direction 
(the director), i.e. no local trihedron of directions, as in the uniaxial nematic A^, the SmA, 
and the columnar D cases." We believe that our work provides the natural extension of the 
theory of continuous defects in solids to the mesomorphic phases named above. 

As a further example of such differences, we note, as mentioned in [KF08J Sec. IIA (p. 
69, Remark), that the Volterra process is properly defined only for magnitude of rotation 
angles less than tt. However, as far as wedge disclinations in liquid crystals are concerned 
the +1 and +| strength defects are entirely different entities - as shown in Section [ij our 
theory is able to predict such differences. 



Computation of bounded energy distortion and di- 
rector solutions 



We demonstrate and evaluate the ideas presented in Section |3.2| through explicit construc- 
tions for different kinds of defects. First, we present results for straight wedge and twist 
disclinations and compare our results with Frank's planar model jFraSSj . Then we take the 
next step forward and simulate twist disinclination loops. 



^Up to a spatially uniform skew- symmetric tensor field. 
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4.1 Straight disclinations 

We refer to disclinations with core cylinder along a straight line, as opposed to a general 
curve or a closed loop, as straight disclinations. A straight wedge disclination is one for 
which the defect line is along the axis of rotation of the director field. In the case of a twist 
disclination, the line is perpendicular to the axis of rotation. 

4.1.1 Wedge disclinations 

Although it is uncommon find an isolated wedge disclination in nature due to its considerable 
energy cost, it is nevertheless a canonical problem for theoretical tools and we study such a 
2-dimensional line defect jKL03 ]. 

Consider as a function of the coordinates Xi and X2, so that the director field would be 



planar. According to (11) and Stokes' theorem 



J (3ldA = - J curl A • CsdA = - ^ X • dr = 2^K. (16) 

This criterion indicates that the line integral of the vector field A over any closed curve 
bounding the core is non- vanishing. Inspired by the Weingarten-Volterra process in elastic 
dislocation theory and adapting it to the simpler situation at hand as indicated in Sections 



3.2 3.4, we define A as a vector field with support on a thin layer that originates from the 



defect core and extends all the way to the boundary: 

X:=J^e, for(f^>0 (17) 

2ri \ \xi\ < ri. ^ ^ 

Here, ri defines the spatial width of the layer. The derivative of A with respect to X2 is zero 
everywhere, except at the layer boundary at X2 = 0. Consequently, will be a singular 
wall-like distribution of finite extent supported on this layer boundary. However, for reasons 



indicated in Section [3^21 the total energy of the disclination so defined would still necessarily 
be finite. 

In agreement with Figure [H we now endow this singular wall with a finite widthj^ Thus, 
the following definition for the layer field is used: 

-^1 \ 2^kr 0+2x2 1^ I <^ro ^ ^1' (^8) 

I 2ri 2ro 1^2 1 ^ 2 

Figure [2] shows the norm of A and the resulting field from it, where both fields have been 
normalized by K: 

A 4a 



^Of course, ideally one would seek such a width as an outcome of the full dynamical model. As a first 
step, we show here that isolated disclinations of specified strength with finite-size cores can be solutions to 
balance of moments in our model. 
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(a) |A| (b) m 

Figure 2: Normalized X field, resulting from (18) and the corresponding field. 



Note that, as mentioned in section 3.2[ the /3 field may not be as locahzed as /3 since 



j3 = cml{E) 


















^ X drti 



and noting that ejkiEfj^ = /3j and ejkiOkOi = 0, 



dd 



(19) 



One can easily see that while the second term is localized inside the defect core, the first 
term is not necessarily zero inside the layer. However, since 6^k decays sharply from the 
defect core, the effect of this term vanishes rapidly away from the core. Figure |3] shows this 
deviation from being a perfectly localized field. 

Equations ([T]), ( [14| ), and (18) yield the director field from the following Poisson's equation 



div {grad 9) = div A. 



(20) 



As a close reference point for the numerically computed approximate director solutions 



that are generated subsequently, the closed-form solution of (20) in an infinite domain is 



derived in the Appendix [A] corresponding to smooth approximating sequences to the discon- 
tinuous layer fields, i.e. ([l7])-(18), we consider. We have used the Galerkin finite element 
method to solve (20). 
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(a) |/3| 



(b) log(|/3|) 



Figure 3: Normalized /3. The affects of the first term in (19) are more obvious in the logarithmic 
plot. The white area in part (h) represents zero for 0. 



In order to compare our energy and director field results with those available in the 
literature (e.g. |KL03j ). we will consider line wedge disclinations in an infinite space. To 
compare the results with Frank's planar model, we consider the following equation with 
Neumann boundary condition: 

div grad 9 = div X onV 
grad 9 • u = grad (f) • u on 

where F is the exterior boundary of the finite domain V and (j) is the angle of the director 
field with the Xi axis in Frank's solution |Fra58j : 

(/) = Ktan-^ +c (21) 

with c a constant. For the purpose of evaluating the boundary condition for the domains 
involved, it suffices to consider grade/) to be given by 

K 

— (-X2ei + xie2) (22) 

(to the neglect of a surface Dirac distribution related to the discontinuity in Frank's solution 
([21])). Moreover, one must assign the value of 9 at one point in the domain for uniqueness 
of solutions. Note that this assignment is tantamount to fixing c in Frank's solution. We 
choose to assign a value ^0 at one point on the right-hand-side boundary of the layer/strip 
on which A has support. 



Now, using this Neumann boundary condition, the weak form of (20) for use in the 
Galerkin finite element method is 

/ 59^,{9^,- X,)dV - [ 59(t)^^Ty,dS = 0. (23) 
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Here 59 is di continuous test function for 9. Bilinear, 4-node, quadrilateral elements are used 
for our simulations. Needless to say, the resulting 9 will be a continuous field and we show 
converged results for the director field with mesh refinement in Figure |15(b)| 

Figure [3] shows the result of solving (23) for 0, plotted as director fields using (|4), for 



different disclination strengths. In all of these examples, ^0 = 1, except in Figure 



4(d) 



These director fields are completely analogous to those from Frank's solution (21), outside 
the layer. 

However, inside the strip, the continuous 9 field has to compensate the jump between 
the two sides of the layer. As a result, one can see that inside the layer, directors rotate 
sharply from one side to the other but at no extra energy cost in the layer beyond what is 
involved in the standard interpretation of the energy due to Frank^s solution, i.e. the energy 



arising from considering the director gradient field to be (22). This is due to compensation 



by the A field - more precisely, by the gradz part of the A field. Moreover, the topological 



content of the standard interpretation, i.e. (22), of the singular gradient of Frank's solution 



(21) is preserved in the E field and approximately in the gradient of continuous 9 field 
(up to the layer). Within this setting, even in the case of an integer defect, although there 
is no apparent discontinuity between the directors on the left and right side of the layer, 
there exists a 27r transition within the layer. In order to have a better understanding of this 
compensation, we take a closer look at the case K = ^. Figure [5 shows the rotation of the 
director field, inside the layer. Note that the rotation of the director field inside the layer is 
opposite in sense to its rotation outside of the field. 

Next, we compare the energy density field of our theory to the corresponding one of 
Frank's model. Using ([5]), the energy term for Frank's planar model will be 



2^2 

(using the standard interpretation of the gradient). Evidently, the total energy in a finite 
body for such an energy density field will not be finite. However, the energy density of the 
defects we model is non singular in most cases and definitely integrable in all cases, so that 
the total energy will be finite. The energy density of our model can be calculated from the 
following expression: 

,jj{n,E) = ^K{e,i-Xi){e,i-Xi). 

Figure [6] shows energy contours normalized with respect to K and outside the defect core: 

In this Figure, we have also compared the decay of the energy density as it moves away 
from the defect core. According to ipp^ such decay must follow a ^ trend. The results from 
our setting shows good agreement with the analytical results outside the core, while being 
nonsingular inside. Figure |6| 

Note that not only is the energy density in our model non-singular at the center, but 
Figure 6(b)| shows that the energy density profile drops dramatically at the center of the 



defect core. Later, in Figure [r4(b)| we will show that the energy level at this point converges 
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(g) 



(h) K = 2 



{i)K 



Figure 4: Computed director fields with wedge disclination at the center for different disclination 
strengths. Part (c) and (d) have the same disclination strength hut different prescribed values Oq. 



to zero. However, this is not a general rule. In fact the energy density profile inside the core 
depends on the distribution of the layer field, A, and consequently /3^. For instance, if we 
cut off the layer field at the center of the defect core instead of letting it decay linearly, i.e. 
use (17) for the layer field definition, the energy density will be as in Figure |7| Figure [7(by 
reveals the fact that although such sharp cut off in the layer field will cause unsymmetric 
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Figure 5: A closer look at the computed director field for K — ^. The colored strip depicts the 
layer field, A. // one makes a clockwise circuit around the defect core, the rotation of the directors 
will he clockwise, outside of the layer, while it is anti- clockwise inside. 






1^0 






lioo 






IJLOO 


A — — Frank's analytical model 








3.2 -0.15 -0 


1 -0.05 


i 0.05 0.1 0.15 0.2 



(a) 



(b) 



Figure 6: (a) Normalized energy density contour outside of the defect core and, (h) The decay of 
energy density when X2 = 0. 



energy density distribution inside the defect core, the energy density wiU stiU be radiaUy 
symmetric outside the core. 

We also note that even though in this case we have a singular wall-like distribution of (3^ 
at the core edge that can potentially cause singularities in the energy density at the sharp 
corners defining the layer, the total energy for a fixed width of layer will always converge as 
long as ri > 0, i.e. the layer has finite width. This is borne out in our computations where 
the total energy corresponding to Figure [7] appears to converge, as shown in Figure |8| 



4.1.2 Straight twist disclination 

Twist disclinations have the axis of rotation of the director field perpendicular to the defect 
line. They are sometimes called perpendicular disclinations. 

From the mathematical point of view, we use the same definition for the planar director 
field as in equation Q. While the n field lies in the Xi — X2 (horizontal) plane, 9 is assumed 
to be a function of Xi and x^ and it is constant along the X2 direction. The definition of the 
layer field undergoes similar, appropriate changes. We define A as in (17), but now in the 63 
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Energy profile on y-axis 
'Energy profile on x-axis 



(a) 



(b) 



Figure 7: (a) Normalized energy density contour when X sharply ends at the center, (h) The 
corresponding energy density profile in cases X2 = (blue, solid line) and xi = (red, dashed line). 




Number of Elements in a fixed, 4 by 4 Domain 



Figure 8: Convergence of the normalized total energy, corresponding to the energy density shown 
in Figure^ with the increase of the number of elements used in a fixed domain. 



direction instead of ei (for direction of A and the width of the layer) as shown in Figure [9(a)| 
We perform a 3-d finite element computation in this case. The computed energy density 
concentration is shown in Figure |9(b)[ The distribution of the energy and its decay on the 
vertical plane is exactly the same as shown in Figure |6] for the straight wedge disclinations, 
which is consistent with expectation. Figures [9(c)| and [9(d) , respectively, show the computed 
and analytical director fields. 



4.2 Twist Disclination loops 

Disclination loops are more commonly found in nature than curves terminating at boundaries 
|HR06j . In this section we will simulate a twist disinclination loop, using our framework. 

The director field will be planar, but now the 6 field will be a function of Xi,X2 and X3. 
The A field for a planar twist disclination loop is taken to be a constant field in a bounded 



layer as can be seen in Figure [10} The bounding cylinder of the layer over which the constant 
field decays to zero then forms the disclination loop. Clearly, the closed loop can be of any 
shape, but restricted to be planar in this exercise. 

In particular, we have chosen a square-shaped and a circular loop in the horizntal plane. 
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0.5 
0.4 




(c) (d) 

Figure 9: Straight twist disclination: (a) 3D plot of the layer field, X, (h) a 3D plot of different 
layers of the normalized energy density to illustrate the energy density distribution in the body, (c) 
is the computed director field, while (d) is the analytical one. 



(3 and (3^ field for the rectangular loop are shown in Figure ]_]_, Just like the 2-dimensional 
wedge disclination case, the (3^ field is non-zero only at the defect core which appears as a 
square loop, and (3 is non-zero inside the loop, where the layer field is non-zero. 

The director field is a planar vector field with twist axis along the direction. As shown 
in Figure [12} in any given plane X3 = constant the effect of the loop decays rapidly with 
distance from the loop. Figure 12 also clarifies the way in which the director vectors make 
a twist as one moves through a circuit enclosing the core. Note that the director vectors 
corresponding to the lower part of the A field have the opposite direction of the ones on 
the top, since that is a half-integer disclination loop (X = |) and as a result the director 
field makes a vr-radian twist. We emphasize that while this twist deficit is observable in the 
(continuous) director field only on circuits that are 'closed up to the thickness of the core', 
integrating on an exactly closed circuit encircling an arm of the disclination loop core 
would yield a 7r-radian twist by design. 
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Figure 10: An example of X field for a square-shaped twist disclination loop. The field vanishes 
in the white region. 





(a) 



(b) 



Figure 11: (a) Normalized for the twist disclination loop and (h) Logarithmic contour plot of 
normalized |/3|. The white area in part (h) represents zero for |/3|. 



Figure [13] shows the energy plots of the two dischnation loops from different points of 
view. Figures 13(a) and 13(b) are the energy density contours at X3 = 0, on the X1 — X2 plane 
for the square-shaped and circular loop, respectively. Note that the energy density of the 
rectangular loop has significantly higher magnitude at the corners, this being a consequence 
of the high curvature in the director field at the corners of the loop. Similar to what we saw 
in Figure 6(b)| these contour plots show a drop in magnitude inside the core. This plot also 
shows smoothness of the energy density, even in the vicinity of the defect core. In order to 



better understand the pattern of the energy density on this plane. Figures 13(c) and 13(d) 



show the logarithmic scale of the plot. A closer look at these two plots reveals that even for 
the rectangular defect core the energy density adopts a radially symmetric patter away from 
the core. The observable deviation from this circular texture on the boundary is due to the 



zero Neumann boundary condition that we have used for these two examples. Figures 13(e] 
and |13(f) show the logarithmic plot of the energy density on the Xi — x^ plane. 
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(a) Square-shaped loop 



(b) Circular loop 



Figure 12: Computed director fields of twist disclination loops for K — ^. 

4.3 Convergence of numerical results 

One of the main advantages of our model is the possibility of attaining non-singular results 
in energy density. In this section we have verified this claim by refining the mesh in our 



numerical model to demonstrate convergence even inside the defect core. Figure 14 shows 
this convergence for the straight wedge disclination. Since the results for diflFerent disclination 
strength are similar, we have normalized our energy plots with respect to K. 
The energy density results for our model and Frank's analytical model are the same, outside 



the core. Figure 15 shows a relative error of less than 1% everywhere for the energy density, 
and less than 0.6% for the director field (outside of the layer) for the half-integer defect. 
Results for all other disclination strengths show the same error distribution. 



Finally, Figure [16] shows the convergence of the energy density in the square-shaped twist 
loop model. In this Figure, the energy density profile of the loop with two different element 
sizes has been compared on the Xi axis, i.e. on the line X2 = 0, X3 = 0. 



5 Balance of forces 

Assuming a single disclination in the body (loop or a core cylinder terminating at the bound- 
aries of the body), static balance of linear momentum ^ is satisfied in the region of the 
body excluding the core. This is so since the stress tensor is given by 

and ([7]) and ^ imply, along with the core-support of /3^, that 

and one chooses the indeterminate constraint pressure field in the excluded-core region of 
the body as 

P= 
20 



(a) 



(b) 




(f) 

Figure 13: Energy plot of the twist disclination for square-shaped and circular loops: (a) and (h) 
are the contours of the normalized energy density, where X3 = 0^ (c) and (d) show the decay of the 
same contours, in a logarithmic scale. Also, (e) and (f) are the logarithmic contour plots of the 
energy density in xi — xs plane, where X2 = 0. 
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Figure 15: Comparison of the computed results from the model with Frank^s model: (a) Relative 
error of energy for wedge disclination, (h) Relative error of director field for wedge disclination. 



to satisfy ([2]). On the external boundary of the body, recognizing (10) one chooses p = 0. 
Thus, the resultant force on the surface of the core cylinder (or any other surface enveloping 
it) must vanish. 

An interesting possibility is afforded by the model |AD12 j in satisfying static balance of 
linear momentum everywhere. Note that within the core curlE^ does not vanish so that it is 
not a priori clear that {E^Ej)j can be expressed as the gradient of a scalar field. However, 
the vanishing of the curl of the field {gradE^)E^ along with the conditions ([7jj8]|9]) within 
the core region may be imposed as a design constraint on the f3^ field. Were solutions to 
exist to this system for /3^, then one could proceed to set up the field A as before with the 
caveat that now the field curlX within the core is not subject only to the mild constraint ([s]) 
but instead to the point wise constraint defined by our calculated field. 



22 




Figure 16: Normalized energy density profile of the twist disclination for the rectangular loop in 
Figure 13(a)\ when X2 — 0. The element sizes are L = 0.02 (red, dashed line) and L = 0.01 (blue, 
solid line). 



Alternatively, the theory also allows the ansatz 

2 r. 

E = 12 Parade- -X-]^ (24) 

and choosing this affords four more degrees of field degrees of freedom to satisfy static balance 
of linear momentum without body forces simultaneously with balance of angular momentum - 
presumably even without necessarily requiring the constraint of incompressibility. Moreover, 
the solutions generated in this paper would seem to suggest that these extra degrees of 
freedom could be invoked only in the core cylinder, thus conferring an automatic multiple 
spatial length scale feature to the theory. Of course, with static balance of linear momentum 
satisfied everywhere (including cores), the resultant force of the Ericksen stress tensor (and 
hence the energy-momentum tensor of the classical nematic theory [E shSOl ICF02j ) on the 
bounding surface enveloping any segment of a core has to vanish without the requirement 
of any externally applied body forces or the onset of fiow, while the |AD12j dissipative 
dynamical theory provides a driving force for disclination motion not related to the Ericksen 
stress or its resultant. This possibility has the potential of reconciling several different 
points of view and clarifying an often confusing debate in the liquid crystal literature (cf. 
|Esh80i lKro93[ lEri95j ) revolving around the physically meaningful question of modeling the 
experimentally observed motion of disclination curves in the presence of negligible or no fiow. 



6 Concluding Remarks 

We have demonstrated the possibility of representing non-singular isolated disclination fields 
within a dynamic theory jAD12 ] of liquid crystalline media designed to account for discli- 
nation and dislocation defects. The solutions satisfy balance of moments everywhere in the 
body, including the core. Balance of forces is satisfied outside the cores and it appears possi- 
ble, by going beyond the ansatz ^ assumed in this work, to satisfy force balance everywhere 
in the body. The fields constructed here may be considered to be equilibria for the theory 
when the defect lines are considered to be externally-pinned (i.e. the defect velocity field, a 
constitutive quantity, is assumed to vanish). 
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Our primary goal has been to demonstrate distinguished pairs of director and (incompat- 
ible) director distortion fields (n, E) satisfying ([T]) and ^ that generate classically familiar 
director and energy density fields of disclinations (but not exactly the same ones) . In doing 
so, we use an intermediary 'layer' field which, thankfully, is not required in the subsequent 
dynamics of the |AD12j model once the director and the director distortion fields have been 
constructed. However, a fact that gives us pause is that the director incompatibility tensor 
fields f3 of these isolated disclinations appear to contain some signature of the layer field 
used to construct them even outside the region that one would typically like to associate 
with the core of the disclination, and whether this is physically realistic or not remains to be 
seen. Also, the director fields that we construct have a rapid variation in the layer (without 
a corresponding energy cost), and the physical realism of such a distribution should also be 
testable in the modeling of situations involving interaction of the director field with fiow. 
Of course, the analysis presented in Section |3.2| provides the logic to construct a director 
and director incompatibility field without the layer-related features discussed above. To wit. 



considering the example of the axial wedge disclination of Section [41] in an infinite domain, 
we define 

Xs = grad Zg + 

where 

div{grad Zg) = div Xg 



and Xg is the distribution obtained in (17) by letting ri 0. This results in Zg being exactly 



Frank's solution (21). Additionally, we let E^ be the solution to 

divE^^ = 
curlE^^ = (3^, 

(with natural decay requirements at large distances) where (3^ is a radially symmetric core 
function described in Appendix with associated solution for E^^. This formula shows 



that outside the core, E^ matches the standard interpretation (22) of the gradient of Frank's 



solution (21). Setting E^ := gradO — Xg and requiring it to satisfy the governing equations of 
our model, i.e. ^ and ([9]), a (discontinuous) solution is certainly 9 = Zg^ with discontinuity 
supported on the half plane defined by ri ^ in the definition (17) of A^. These results 



imply that outside the core and away from the half-plane Xi = 0, X2 > 0, E^ matches grad 9. 



Assuming the validity of (15) on the half-plane and noting that the grad 9 and E fields are 
parallel to ei on it, it may be concluded that /3 is non- vanishing only inside the core. With 
E^ defined as in this paragraph, we note that even though grad 9 is not a locally integrable 
field, this is exactly compensated by gradzg so that the energy density is formed from the 
field E^ alone. In this connection, we note the interesting fact that the classical defect theory 
works with a director- 'gradient' field for forces, moments and energetics that is not strictly 
a generalized derivative of what it states as the director field (even outside the core 

fl Thus, 

^The analog of Weingarten's construction to be utilized here makes it entirely clear why this must be so - 
there, one begins with an irrotational field in a doubly-connected 3-d domain with a through-hole and what 
may be associated as the 'closest' potential field generating the irrotational field is not defined on the same 
domain but rather on a simply-connected domain obtained through a single cut of the doubly-connected 
domain that renders the latter simply-connected. Unfortunately, however, these are not procedures that 
lend themselves to practical nonlinear theory, especially of the type interpretable by the intelligence level of 
electronic computers that are essential to probe the complex dynamics of the theories we have in mind. 
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it would seem that considering the energy formed from a director distortion field distinct 
from the director gradient in the presence of line defects (as we do) is probably a reasonable 
logical device even for discussing the classical theory. 

This paper has not made any statements about whether the chosen isolated defect-like 
director and director distortion fields can actually occur as equilibria of the overall dynamical 
model |AD12j for suitable constitutive assumptions for defect velocities (e.g., a linear kinetic 
assumption on disclination velocity as a function of its theoretical driving force). Based on 
our experience with related work on Field Dislocation Mechanics |AchlQ| , |AT11| , jAMZlO] , 
we believe this should be possible with the use of non-convexity in OZF energy (possibly 
due to a limit on director elasticity), the system (as opposed to scalar) nature of the defect 
evolution equations, as well as from appropriate boundary conditions. In the presence of 
fiow, there are, of course, other natural means to induce line defects. 

There are many natural further steps to explore: solving for more-than-planar fields of 
isolated defects, defect solutions for the full 4-constant Oseen- Frank energy, solutions for 
disclinations coupled to dislocations in smectics and columnar phases and on to dynamics 
with and without fiow (without fiow and even without director inertia, the model has dynam- 
ics due to defect evolution). While each of these problems is significant in its own right, we 
take some satisfaction in the fact that all of these questions can at least be posed unambigu- 
ously within our modeling framework. Finally, an important future task is to make contact 
with the fine topological work of Kleman and Friedel |KF08| and Kamien and co-workers 
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Appendices 

A An Explicit solution for 6 



An analytic solution for the director field corresponding to A source fields in [20] in an infinite 
domain is developed, to serve as motivation for some of our claims in Section [3^ The source 
fields are layer-like but do not have compact support on bounded domains. Consider a 5 
sequence |Han06j : 

5n{x) 



So that, j^^Sn{x) = erf(x) = 1, for any n. Also, consider a H sequence, which converges 
to the Heaviside step function 

Hn{x) = - (tan~'^(nx) + l) . 
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Now, one can define A as 

Xn{x) = 27lK5n{xi)Hn{x2)ei. 

Note that, for any given n. 



j) Xn.dr = 27rK J 5n{xi)dxi = 2tiK 



In order to solve (20), one can use Green's function method for an infinite 2D domain, where 
the Green's function is 

G{x,x') = ^H\x-x'\). 
Now, the expression for 6^ in terms of the Green's function wiU be 

9^= I G{x,x')dwXn{x')dV 
Jv 

Integrating by part and considering the fact that the field has only the ei component 
9^= ! G{x, x')K{x') -udS- [ -j^G{x, x')K{x')dV 



Jyf \x' — X\ 



./a:; Jx'^ |x' - X 



27rk5n{x[)Hn{x2)dx2dx[ 



As a result: 



9^ = ^ [ [ ^^^e-(^^i)' (tan-i(nx'2) + l) ^^^^i (25) 
2v^ Jx[ Jx'^ \x' -x\ 



Note that in case of n ^ oc. 



9oo = ^ [ [ — ^5oo{x[)H{x2)dx2dx[ 
^ Jx' Jx' \x' - x\ 



K r ^1 



2 Jx'^=0 ^1 + (^2 X2 
K _T fX2 — X'<2 

= — tan ' 



2 V xi / 4=0 

as a result, 9n will converge to 

^ = Xtan-^ +c, (26) 

which is exactly Frank's solution for planar director fields for straight disclinations for the 
1-constant energy. 
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e 



B Solution for E, 

We provide the following solution from |AD12j for completeness. Consider /3f given by 



with the stipulation that 



0,r > Tr 



(27) 



2K7T 



277 rr, 



Jo 



p{r) dr rdijj ^ K 



p{r) rdr 



(28) 



where K is the strength of the disclination and Vc is a core radius. Then following Section 
5.2 in [ AchOlj . we have the solution 



For r < Tc : E^^ 



-X2 



p{s)s ds, E% 



Xi 



p{s)s ds 



For r>r.: ^f, = i^^, E',, = K^. 



(29) 



It can be checked that this solution indeed satisfies the governing equations for for any 
p satisfying the given conditions. For the choice 



p{r) = 

the inside-core solution evaluates to 

. -X2 f2K 
For r<r,: Ki = -^[ — 



2K f 1 1 

0, r > 



2rr 



xi f2K 



2rr 



while the outside-core distribution remains unchanged. 
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